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Abstract 

Trapped ions are considered one ol the best candidates to perform quantum information pro- 
cessing. By interacting them with laser beams they are, somehow, easy to manipulate, which makes 
them an excellent choice for the production of nonclassical states of their vibrational motion, the 
reconstruction of quasiprobability distribution functions, the production of quantum gates, etc. 
However, most of these effects have been produced in the so-called low intensity regime, this is, 
when the Rabi frequency (laser intensity) is much smaller than the trap frequency. Because of 
the possibility to produce faster quantum gates in other regimes it is of importance to study this 
system in a more complete manner, which is the motivation for this contribution. We start by 
studying the way ions are trapped in Paul traps in and review the basic mechanisms of trapping. 
Then we show how the problem may be completely solved for trapping states; i.e., we find (exact) 
eigenstates of the full Hamiltonian. We show how in the low intensity regime Jaynes-Cummings 
and anti-Jaynes-Cummings interactions may be obtained, without using the rotating wave ap- 
proximation and analyze the medium and high intensity regimes were dispersive Hamiltonians are 
produced. The traditional approach (low intensity regime) is also studied and used for the gener- 
ation of non-classical states of the vibrational of the vibrational wavefunction. In particular, we 
show how to add and subtract vibrational quanta to an initial state, how to produce specific super- 
positions of number states and how to generate NOON states for the two-dimensional vibration of 
the ion. It is also shown how squeezing may be measured. The time dependent problem is studied 
by using Lewis-Ermakov methods, we give a solution to the problem when the time dependence 
of the trap is considered and also analyze an specific (artificial) time dependence that produces 
squeezing of the initial vibrational wave function. A way to mimic the ion-laser interaction via 
classical optics is also introduced. 
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1 Introduction 



The possibility to trap small clouds of particles, or inclusive to trap individual atoms or ions, in an 
small region of space, was opened with the invention of electromagnetic traps. These traps allow to 
study isolated particles for long periods of time. The Kingdon trap [1] is considered the first type of 
trap developed, it consists of a metallic filament surrounded by a metallic cylinder, and a direct current 
voltage applied between them; the ions are attracted by the filament, but its angular momentum makes 
them turn around in circular orbits, with a low probability to crash against it. A dynamic version of 
this trap can be obtained if also an alternate current voltage is applied between the poles. However, 
this type of trap was not widely used at that time, because it has short storage times and because its 
potential is not harmonic. In 1936, Penning invented another trap [2] in which the action of magnetic 
fields together with electric fields make possible the trapping of ions. The complete development of 
this type of trap was reached when, in 1959, Wolfgang Paul designed an electrodynamic trap (now 
called Paul trap) [3]. In the Paul trap the idea is that a charged particle can not be confined in 
a region of space by constant electric fields, instead an electric field oscillating at radio frequency, 
must be applied. The Paul trap uses not only the focusing or defocusing forces of the quadrupolar 
electric field acting on the ions, but also takes advantage of the stability properties of the equations of 
motion. The ions trapped individually are very interesting, mainly because they are simple systems 
to be studied. In particular, we take advantage that the ion motion in the Paul trap is approximately 
harmonic, making this system a simple one, allowing a better and more direct comparison with theory. 
Individual ions of Ca+, Be+, Ba+ and Mg+, can be storage, even for several days. The trapped ions 
can be used to implement quantum gates, and a bunch of ions arranged in a chain, is a promising tool 
to achieve a quantum computer (each ion in the chain is a fundamental unit of information or qubit) 
[4] . The trapping of individual ions also offers a lot of possibilities in spectroscopy [5J , in the research of 
frequency standards [SI [7] , in the study of quantum jumps [5] , the engineering of specific Hamiltonians 
[9] and in the generation of nonclassical vibrational states of the ion [TOl [iTJ [121 [ISl [HI [HI HB [O] , to 
name some. To make the ions more stable in the trap, increasing the time of confinement, and also 
to avoid undesirable random motions, it is required that the ion be in its vibrational ground state. 
This can be accomplished by means of an adequate use of lasers; with the help of these lasers, the 
internal energy levels of the trapped ion can be coupled to their vibrational quantum states, in such 
a way, that for a certain detuning, the coupling is equivalent to the Jaynes-Cummings Hamiltonian 
[T8l fT9l [20l [2TJ [22l [23]. On the other hand, the beam that induces the coupling can be tuned to 
allow interactions that generate simultaneous transitions of the internal and vibrational states, either 
to lower vibrational energy levels (while passing from the excited to the ground state) or to higher 
vibrational energy levels (while passing from the ground to the excited state); this type of coupling is 
called anti- Jaynes-Cummings. Alternating successively Jaynes-Cummings and anti- Jaynes-Cummings 
interactions, the trapped ion can be driven to its vibrational ground state. In this report, we study part 
of the physics of the trapped ions interacting with a laser field. By using a set of time-dependent unitary 
transformations, it is shown that this system is equivalent to the interaction between a quantized field 
and a two level system with time dependent parameters. The Hamiltonian is linearized in such a way 
that it can be solved with methods that are found in the literature, and that involve time-dependent 
parameters. The linearization is free of approximations and assumptions on the parameters of the 
system as are, for instance, the Lamb-Dicke parameter, the time-dependency of the trap frequency 
and the detuning; thus, we can obtain the best solution for this kind of systems. Also, we analyze a 
particular case of time-dependency of the trap frequency. 

Because we will assume an ion trapped in a Paul trap, in Section 2, we review the basic mechanisms of 
it. In Section 3, we show how exact eigenstates may be obtained for the ion laser system. In Section 4, 
we show how this interaction may be solved in different regimes. In Section 5, it is treated the standard 
approach to this interaction; i.e., the low intensity regime, where by application of the rotating wave 
approximation, Jaynes-Cummings and anti- Jaynes-Cummings interactions may be generated. By using 
this approach, we show how nonclassical states may be generated; in particular, we show how to add 
phonons to a vibrational state, how to filter specific superpositions of the motional wave function and 
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how to generate NOON states in ions vibrating in two dimensions. We also propose here a scheme that 
allows the measurement of squeezing. In section 6, we analyze the case of an ion with a time-dependent 
frequency interacting with a laser beam. By doing a series of unitary transformations, we linearize the 
Hamiltonian of the system to an exact soluble form; this linearization is also valid for any detuning 
and for any time dependence of the trap. In Section 7, we show how the ion-laser interaction may be 
modeled by evanescently coupling waveguides and Section 8 is left for conclusions. 

2 Paul trap 

2.1 The quadrupolar potential of the trap 

As we already said, the Paul trap uses static and oscillating electric potentials to confine charged 
particles. A charged particle is linked to an axis if a linear restoring force acts over it; i.e., if the force 
is 

F = -cf, (1) 

where f is the particle position and c is a constant. In other words, if the particle moves under the 
action of a parabolic potential, that can be written in general form as 

<^{x,y,z)^A{ax^ +I3y^ +jz^), (2) 

where A is another constant. The potential $ must satisfy the Laplace equation, which means that 

V2$==0, (3) 

where is the Laplacian operator. The Laplace equation ([3]) imposes the condition 

a + /3 + 7 = 0. (4) 

To satisfy the above condition, we have several possibilities: 

a) We make a = 1, /3 = and 7 = — 1, and this takes us to the bidimensional potential 



b) Another possibility is a — 1, 13 — 1 and 7 — —2, and in this case we have, in cylindrical coordinates, 
the potential 



'^ = ^ai^ir^~^-')^ (6) 
'0 + ^^0 



with Tq = 2zg. 

Configuration a). Equation ([S]), is created by four hyperbolic electrodes linearly extended in the z 
direction, as shown in Figure 1. Configuration b), expression ([5]), is created by two electrodes with the 
form of an hyperboloid of revolution around z axis. 

The most used trap is the linear trap, as the one shown in Figure 1, but with poles with circular 
transverse section instead of hyperbolic, because it is easier to build. This cylindric form does not 
correspond to some set of values of ([1]), but numerically it has been demonstrated that the potential 
produced by these electrodes near the axis of the trap is very similar to the one produced by the 
hyperbohc electrodes [24] . 

For the tridimensional case the magnitude of the field is given by 

Ex = —X, Ey = —y, E^ = 2—z. (7) 







Expressions ^ and ([7]) reveal that the components r and z of the electric field are independent from 
each other, and that they are linear functions of r and z, respectively. We also see that we have a 
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Figure 1: Electrode structure for the bidimensional configuration expressed by Equation ([S]). 

harmonic oscillator potential (parabolic and attractive) in the radial direction and a parabolic repulsive 
potential in the z direction. If a constant voltage $0 is applied, and an ion is injected, the ion will 
oscillate harmonically in the x — y plane, but because the opposite sign in the field Ez^ its amplitude 
in the z direction will grow exponentially. The particles will be out of focus, and they will be lost 
by crashing against the electrodes. Thus, the quadrupolar static potential, by itself, is not capable to 
confine the particles in three dimensions; at most, with this potential, we get unstable equilibrium. 
We will see next, how to solve this problem. 

2.2 Oscillating potential of the trap 

To avoid the unstable behavior of the charged particles under a static potential, the trap must be 
modified. If an oscillatory electric field is applied, the particles can be confined. Because the periodic 
change of the sign of the electric force, we get focusing and defocusing in both directions of r and z 
alternatively with time. 

If the applied voltage is given by a continuous voltage plus a voltage with a frequency SI, we have 



'^o = Uo + Vq cos Vlt, 



(8) 



and the potential in the axis of the trap is 



$ = 



Uq + Vq cos Q.t 



{r'~2z^), 



(9) 




where rp is the distance from the trap center to the electrode surface. 
In Figure 2, we show a transversal section of a Paul trap using an oscillating electric field. 



4 




e 



7 



r ® 




IL 



(\ 



Figure 2: Scheme of a Paul trap to storage charged particles using oscillating electric fields generated 
by a quadrupole. The figure shows two states during an alternate current cycle. 



2.3 Motion in the Paul trap 

We will study now some details of the ion motion in a Paul trap. Let us consider the particular case of 
just one ion, in three dimensions. If m is the mass of the ion, and e its charge, the equation of motion 
is 

mf{x,y,z) = qE = -qV^. (10) 
In order to analyze the trapping conditions, we write explicitly each component, 

2e 



2e 



and 



where 



2e 



{Uo + Vo cos fit) X, 
{Uo + Vocosnt)y, 



mB? 



2(C/o + yocosm)z, 



2zl 



Making the substitution 



8eC/o 



' 2 



qr 



46^0 



' 2 



nt 



equations ([TT|) to (fO|) . take the form of the Mathieu equation; i.e., take the form 

X 

+ (ar — 2qr cos 2t)x — 0, 



and 



cPy 

— r + [ar - 2qr COS 2T)y = 0, 



—J + (a^ - 2qz COS 2t)z = 0, 



respectively. We can write the three equations as the following one, 

(fu. 



+ {ai — 2qi cos 2r)wi = 0. 



(11) 
(12) 

(13) 

(14) 

(15) 
(16) 

(17) 
(18) 



The subindex i — r, z corresponds to the quantities associated with the axial and radial motions of the 
ion, respectively. The quantities Ui represent the displacement in the directions r and z. 



5 



2.4 Approximated solution to the Mathieu equation 

The Mathieu equation is a hnear ordinary differential equation with periodic coefficients. This equation 
can be solved using Floquet's theorem [25 , which takes us to the general solution 

u,{t) = A,e"5-^0(r) + B,e-'^'^(/.(-r), (f9) 

where Ai , Bi and jSi are constants determined by the initial position, by the initial velocity of the ion, 
and by the trap parameters a and q, and 

+00 

0(T)=0(r + ^)= ^"^2"" (20) 

n— — 00 

is a periodic function. 

The Mathieu equation has two types of solutions: 

f) Stable motion. When the characteristic exponent (3 is real, the variable u(r) is bounded, and in 
consequence the motion is stable. That means that the particle oscillates with bounded amplitudes 
and without crashing against the electrodes. These conditions allow to trap the ion. 
2) Unstable motion. When the characteristic exponent /? has an imaginary part, the function u{t) has 
an exponential growing contribution. The amplitudes grow exponentially and the particles are lost 
when they crash against the electrodes. The boundaries of the stability regions correspond to zero and 
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integers values of Pi, and the first region of stability is surrounded by the four lines Pr — ^ , Pr ~ 
Pz = 0, and /3z = 1, as is shown in Figure 3 

As Pi is determined by a and q, the Mathieu equation has stable solutions as a function of a and q. 
Stability regions for the solutions of Equation ()18|) correspond to regions in the space of the parameters 
a — q, where there is an overlap of the stability regions in the axial and radial directions. 
In the literature it is not possible to find analytic solutions for Equation ()18|) . but in most of the 
applications an specification of the map of stability of the solutions is enough, and it is not necessary 
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a detailed functional dependence. However, an approximated solution can be given in the stability 
region of interest. To this end, we can write expression ([T^ as 



+00 +00 



u,{t)^A, ^ C^„cos(2n + A)r + B, ^ q„sin(2n + A)T (21) 

n— — oo n— — 00 

where, as we already said, Ai and Bi are determined by the initial position Ui{Q) and initial speed 
iii (0) of the ion, respectively. The subindexes i = r, z coincide with the quantities associated with the 
radial and axial ion motion, respectively. 

The coefficients in the solution ([2T|) are given by the recurrence relations 



^2n+2 ~ ^2n^2n + ^'2,1-2 — (22) 

with 

ni„ ^ (23) 

Once given and qi, the quantities C|„ and /3i can be calculated. If we define, 

(Ji 

^2.1 = ^i^ACQ, B[^BiCQ, (24) 

and we make 

u.[t) = ul{t)+uT{tl (25) 

we get, from Equation ([2T|) . 

(r) = A ■ cos + B[ sin w^i , (26) 

and 



u™(t) = ^(A-coswit + S-sina;jt)(G''^„ + Gian) cos nfJi 



(27) 



+ (i?-cosi:jii - A- sin Wit) (Can - G"_2„) sinn^it. 



where Ui = l3iQ/2. 

Analyzing Equations (|26p and ([77]) . it is possible to realize that the ion motion has two components: 
uf{t), a harmonic oscillation of frequency LOi, and, u^{t), a superposition of several harmonics with a 
fundamental frequency 17, and amplitudes modulated by the frequency uji. However, the proportion of 
the two components, the values of Wi, the number of subcomponents that contribute appreciably and 
their weights, depend strongly on the values of ai and qi, in such a way that they will change in the 
stability region. All parameters are determined when the values /3i and G^n are given. Several values 
of /3i and G2„, corresponding to some typical values de Ui and q^, are listed in [27] . In the table there, 
it is possible to see that in the first region for a ^ q <C 1, we can assume that G2 = G!_2 and the rest 
of the coefficients G!|_2„, n > I, can be ignored; thus. Equations ([^5]) and (P7)) can be rewritten as 

(r) = Uio cos{ujit + 6i), (28) 

and 

u™(r) = CuiQ cos fit cos (wji + Jj); (29) 

or in other words, 

Ui{t) = UiQ cos{ijJit + Si) {1 + C cos V,t), (30) 

with 

w,o^=Af+i?f, (31) 
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Figure 4: Micro motion and secular motion of a trapped ion with parameters q = 0.2, /3 — 0.02. The 
oscillations at high frequency are the micro motion and those at low frequency are the secular motion. 




5 ^ arccos ' , (32) 



and C = 2Gi. 

In Figure 4, we plot Equation ([30]). The ion motion is composed by two types of oscillations: the 
harmonic oscillation with frequency uii, called secular motion, and the small contributions oscillating 
at the frequency fi, called micromotion. Usually, the micromotion is ignored, but it can be reduced 
using additional electrodes [28j . In this way, the ion motion is controlled by Equation and behaves 
as it was confined in a harmonic pseudo-potential, that for the radial part, has the form 

$, = ^(0.^x^+0;^). (33) 

Typically Uq — 0, thus a = (Equation (IT4)) ): in any case, we are working in the region where a ~ 0. 
Thus, the frequencies uix and iVy are degenerated, and Equation (|33|) reduces to 

= ^ (x' + y') . (34) 
To obtain an expression for w^, we can use the approximation |29) . 



f3r^\lar + f, (35) 



with the definition cj^ = /3rf^/2, to get 



23/2 ^^r^Q- 



(36) 



Experimentally, the typical ranges of operation are Vq ~ 300 — SOOVolts, r2/27r w 16 — ISMHz, and 
ro ~ 1.2mm, that gives a radial frequency ojr ~ 1-4 — 2.0MIIz for calcium ions C'^Ca^). 
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We can summarize this section, saying that under certain conditions it is a good approximation to treat 
the ion motion as a harmonic oscillator. This allows us in the following to apply algebraic techniques 
well known for this system to the ion-laser interaction. In particular, invariant methods |30[ I31i 132) 
that permit the handling of time dependent systems in a simple way. 

3 Ion-laser interaction in a trap with time-independent fre- 
quency 

Despite the relative simplicity of the ion-laser interaction, the full theoretical treatment of its dynamics 
is a nontrivial problem, because this interaction is highly nonlinear. Even in the simplest case, where 
only a single ion is in the trap, one is usually forced to employ physically motivated approximations in 
order to find a solution. A well-known example is the Lamb-Dicke approximation, in which the ion is 
considered to be confined within a region much smaller than the laser wavelength. Many treatments 
also assume a weak coupling approximation; i.e., a sufficiently weak laser- ion coupling constant. Under 
these conditions, tuning the laser frequency to integer multiples of the trap frequency results in effective 
Hamiltonians of the Jaynes-Cummings type jlSj |19j [20l [211 EH ES], in which the center-of-mass of the 
trapped ion plays the role of the field mode in cavity QED. Recently, a new approach to this problem 
has been suggested, based on the application of a unitary transformation T that linearizes the total 
ion-laser Hamiltonian. Under this transformation the Hamiltonian becomes exactly equivalent to the 
classic Jaynes-Cummings model (JCM), including counter- rotating terms and an extra atomic driving 
term. Remarkably, the ion-trap system is thus formally equivalent to an atom interacting with a 
single-mode quantized electromagnetic field. 

This equivalence does not of course lead to an exact solution to the ion-trap problem, since the 
complete JCM is also notoriously unsolved. To our knowledge, eigenstates and eigenvalues for this 
model have only ever been obtained either numerically or expressed in the form of series with no 
known closed expression. Still, some advantage can be taken of regimes where the JCM is analytically, 
albeit approximately, soluble, such as the well-known "weak-coupling" limit where the rotating-wave 
approximation can be made. Using the already mentioned map to translate this solution back into 
the ion-trap scenario has led to the identification of a new soluble regime for that system. This in turn 
has already proven useful in the context of quantum computing. 

In this section, we study the interaction of a laser with a trapped ion in a harmonic potential with 
constant frequency. We start with the Hamiltonian of the system, and we show that it is possible 
to find Jaynes-Cummings type transitions and anti- Jaynes-Cummings type transitions, depending on 
the different cases of resonance and laser intensities that induce the coupling between the ion internal 
states and the ion vibrational states. 
We can write the Hamiltonian of the trapped ion as 



where H^ih is the ion's center of mass vibrational energy, iJat is the ion internal energy, and iJjnt is 
the interaction energy between the ion and the laser. As we explained in the previous section, the 
vibrational motion can be fairly approximated by a harmonic oscillator. Internally, the ion will be 
modeled by a two level system. In the interaction between the ion and the laser, we will make the 
dipolar approximation, so we will write the interaction energy as — er • E, where — er is the dipolar 
momentum of the ion and E is the electric field of the laser, that will be considered a plane wave. 
Thus, we write the Hamiltonian explicitly as 



H — Hvih + -ffat + Hi 



(37) 




(38) 
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The first term in the Hamiltonian is the ion vibrational energy; in the ion vibrational energy, the 
operator n = a^ais the number operator, and the ladder operators d and d^ are given by the expressions 



and 




where we have made the ion mass equal to 1. Also, for simplicity, we have displaced the vibrational 
Hamiltonian by 1^/2, the vacuum energy, that in this case is not important. 

The second term in the Hamiltonian corresponds to the ion internal energy; the matrices cr^, 17+ , and 
c7_ are the Pauli matrices, and obey the commutation relations 

[az,(J±] = ±2a±, [ct+,o-_] = 0-2, (41) 

and LO21 is the transition frequency between the ground state and the excited state of the ion. 
Finally, the third term, is the interaction energy between the ion and the laser; in this last term, we 
have used again the rotating wave approximation. 



3.1 Exact eigenstates 

In this subsection we show that, under certain combinations of system parameters, it is possible to 
obtain exact eigenstates for the ion-trap Hamiltonian. Using the map T we also obtain therefore exact 
eigenstates of the complete JCM. The set of states we construct is by no means complete, and it is 
also unclear at present how (or if) it may be extended. Nevertheless, we believe that its existence may 
provide a clue to a deeper understanding of both models. 

Let us start by recalling the equivalence between the ion-trap system and the JCM |18 [fT9l[20l[2n[22l 
1231 1331 . The Hamiltonian for the ion-laser interaction can be written as 



[a+D{trj)+a-D\tr^)), (42) 



where D{ia) = e*"(''+°*) [34| is the displacement operator, v is the harmonic trapping frequency, 
S — UJ21 ~ the laser-ion detuning, the (real) Rabi frequency of the ion-laser coupling and 



rj = Kd-^ (43) 
V Zmh' 

is the so called Lamb-Dicke parameter, that is a measure of the amplitude of the oscillations of the 

ion with respect to the wavelength of the laser field represented by its wave vector K. 

On the other hand, the Jaynes-Cummings Hamiltonian with counter-rotating terms is given by 

HjcM = + (T^ 4- iA ((T+ -f (T_ ) (a - a^) . (44) 

Although these two models appear to be physically and mathematically quite distinct, they are in 
fact exactly equivalent. The easiest way to see this is by rewriting Equation (|42l) in a notation where 
operators acting on the internal ionic levels are represented explicitly in terms of their matrix elements, 
as 



z/n -I- I flD [ir]) 
QD^ [ir]) 



(45) 



Consider now the unitary operator 

"^^ 7! V -D'(/3) b{fi) 



10 



where (3 = irj/2. It is possible to check after some algebra that 

nj. _rr.rr. rj.^ _ ( vfl + ^ + ^ 
I Lion — ion — I / r \ c- 

\ ^ (a - at) + I 
Returning to the usual notation, we obtain 

•Hion = vn + na,+ '^ (cr+ + cr_) (a - o)) + {a+ + + (48) 

Comparing with Equation (|44| . it can be seen that this is precisely the Jaynes-Cummings interaction, 
supplemented by two additional terms: the first corresponds to an extra static electric field interacting 
with the atomic dipole, and the second is just a constant energy shift which can be disregarded. In 
particular, a purely Jaynes-Cummings form is recovered when 5 — Q, corresponding to a resonant laser- 
ion interaction in Equation (j42l) . Of course, in Equation (|48p the various parameters in the Hamiltonian 
have different meanings than they do in Equation (|42l) : u becomes the cavity field frequency w, 251 the 
atomic transition frequency wq- <^ the coupling strength with the static field and -q the ratio between 
the Jaynes-Cummings Rabi frequency 2A and the cavity frequency uj. In what follows, we shall refer 
to Equation ([48]) as the 'Jaynes-Cummings picture' of the ion-trap Hamiltonian, Equation (|42|) . 
This correspondence is very useful, since it enables one to map interesting properties of each model 
onto their counterparts in the other. For instance, it has been recently used to identify the existence 
of "super- revivals" in the ion-laser interaction [33j , and to discover a means of realizing substantially 
faster logic gates for quantum information processing in a linear ion chain [35| . In this paper, we will 
use it to map eigenstates of one model into those of the other (it is clear that, if ji/') is an eigenstate 
of ifion, then T is a corresponding one for Hion)- In this regard it is important to point out 
that, although Hion and 'Hion, both describe systems consisting of a two- level atom interacting with 
a bosonic mode, one should not identify each of these subsystems with their counterparts after the 
transformation has been applied. This is due to the fact that T is an entangling transformation: ion- 
trap states, where the ion's internal and vibrational degrees of freedom have well-defined pure states, 
can be mapped into entangled atom-cavity states in the corresponding cavity QED system. 



(47) 



3.1.1 Simple ansatz 

Let us return now to the ion-trap Hamiltonian, Equation (I42[) . We will construct an ansatz which 
allows the determination of exact eigenstates of this system, provided certain relations are satisfied 
between the parameters fi, S, and rj. In order to motivate our general solution, let us consider first the 
possibility of finding such a state of the form 

|^) = |e)(co|O)+ci|l)) + |.9)|0)^ f ^"|0) +ci|l> ^ (49) 



where again we choose a notation where the ionic elements are written out explicitly (e.g., |e) = (o))- 
Let us now see whether the eigenvalue equation 

Hionl^) = E\^b), (50) 

can be satisfied. Equation (|^5|) shows that this requires |0) to be of the form 

10) = I?t (iry) (do |0) + di |1» = do \-iri) + di |-i?7; 1) , (51) 
where l—iij) is a coherent state and \a;k) = D (a) \ k) is a displaced number state [36]. We thus require 

( (cof + Mo)|0) + (Mi+ci(z.+ f))|l) \ 
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Now, using the well-known fact that {a) aD {a) = a + a [34], it is easy to show that displaced 
number states satisfy the recursion relation 

h\a]k) = (|a|Vfc) \a]k) + a%/fcTT |a; /c + 1) + |a; fc - 1) . (53) 

Substituting then Equations (|49)) . (|52l) and (l53l) into Equation ([50)) gives the following eigenstate 
conditions: 

VL iriv 5 

di=0; co = -; ci = E = ly + (54) 

which hold however only if the parameters ^2, i5, 77 satisfy the additional constraint 

- +772 = 1 + -. 55 
v J V 

Under these conditions the state 

k+„> = V) |0) + ^ |i)) + 1.9) ■ (56) 

is an (unnormalized) eigenstate of -ffion with eigenvalue v + (5/2. Using operator T we can map this 
state into an eigenstate of the generalized JCM model in Equation 



kjcM> = T |7/^+„) = I-) + ^ 1)) + |+) . {57) 

where|±) = (|5)±|e))/V2. 

Condition (j55p means that the ansatz in Equation (|49p does not always succeed, as only two of the three 
parameters il,,S,ri can be chosen independently. In addition, the domain of some of these parameters 
is not entirely unrestricted; for instance, it is easily seen that no solution exists when the laser is 
tuned such that 6 < -~2v. Nevertheless, the existence of solutions satisfying Equation ([49| leads us 
naturally to seek for other solutions using similar or slightly generalized ansatz. For example, a second 
solution can be easily found if we note that the Hamiltonian i/ion is invariant under the combined 
transformations 

|e) o (5 O 770-77. (58) 

Applying this symmetry transformation also to Equations ([54]) and ([55l) we can see that, as long as 
we satisfy the condition 

+ 7^2 = 1--, (59) 



then 



|V'i;„) = |e)|777) + |.g)(^|0)-^|l)) (60) 



is an eigenstate with eigenvalue u — 5/2. 

Note that, unless 5 — Q, conditions (|55p and (1591) are mutually exclusive. (In the 'ion-trap' picture, 
this means the laser must be resonant with the ion; in the 'JCM' picture, that no static field is applied; 
i.e., that the Hamiltonian is just the full JCM). In other words, only in this special case are iTAion) 
and jV'ion) simultaneous (in fact, degenerate) eigenstates of -ffion- The reason is that, in this case 
only, the ion-trap Hamiltonian, Equation ([42]), has an extra symmetry: it commutes with the operator 
(Ta: exp(z7ra^a). This parity-like observable has two eigenvalues ±1, and so the spectrum of i/ion is 
two-fold degenerate (in the JCM picture, the corresponding symmetry operator is exp(i7ra^a)) . It 
can be easily checked that neither jV'k^n) '^'^'^ |V'ion) have this symmetry, but simple linear combinations 
of them do. 



12 



3.1.2 General ansatz 



One can easily generalize Equation (l56t to obtain a more general eigenstate for Hjon', this general state 
can be proposed as 

^ m+l m 

1^) = - E c» I") 1^) + E I-"?' Is) • (61) 

n=0 n=0 

Substituting this proposal in the eigenvalue Equation (j50p . it can be shown that the corresponding 



eigenvalues are E„i = (m + + |, that 



m+l— n 

2 

1 



n = 0, 1,2, m 
^; n = m + 1 



(62) 



and that the c?„ coefhcients satisfy 



£0 -177 

irj El —i\/2rj 



i\/2r7 



where 



do 

dm 

2 ^ 1 fi2 

e„ = l + m- n- 77 H — . 

V 1 + m — n ly"^ 



iy/mi] Era 



(63) 



Equation (|63p establish the recurrence relations 

iVjildj-i + Sjd-j - iy^j + li]dj+i =0, j = 0, 1, 2, m - 1 



(64) 



(65) 



that allow us to calculate all the d coefficients of the eigenvector (|6T|) in terms of do. The first ones 
are given by 



di 



-leo- 



do 



(66) 



, / 2 \ do 

d2 = iv - eoe^) ^ 



di = (377"* - 3?7^eoei - 2T]'^eoe3 - 77^£2e3 + eo£ie2e3) 



do 



Note that the vector of coefficients (do, ...,dm) is an eigenvector of this tridiagonal matrix with zero 
eigenvalue. This is only possible if the determinant of the matrix in Equation (|63p is zero, which 
imposes a constraint on il, S, v. These imposed conditions are the generalization of Equation (1551) . 
It is easy to show that if £7 ^ t^, ry <C 1 and 5 = —v all these conditions are satisfied (the matrix in 
Equation (j63p becomes diagonal with a zero in the m + l,m + 1 position), and then we obtain the 
exact eigenvalues and the exact eigenvectors without using the rotating wave approximation. 
As we already mention, the Hamiltonian is symmetric under the combined transformations 
{|e) O Ig) ; (5 o ~8\ 7] o —rj}; this allow us to propose another set of eigenstates as 



n 



m+l 



n=0 



dn n) |e) + - ^ c„ \n) \g) , 



(67) 
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where now 



-1- 



-irj^y/rn + Tdfn] n ~ m + I 
The d coefficients satisfy Equation (|63|) but witli 

s 1 



0,1, 2, ...,m , , 

(68) 



V 1 



V\ (69) 



and the corresponding eigenvalue are = (m + l)j^ — |. 

In this case all the constraints imposed by Equation (1551) are satisfied if ^ z^, 77 <C 1 and (5 = and 
the exact eigenvalues and the exact eigenvectors are again obtained without using the rotating wave 
approximation. 



3.2 Blue and red sidebands 

Processes where simultaneously the internal excitation and the motional quantum numbers are in- 
creased (decreased) are known as blue sideband excitations [37]. While, when simultaneously the 
internal excitation is excited (lowered) and the motional quantum numbers are decreased (increased) 
are known as the red sideband. We discuss now some properties of the eigenstates ((6T|) and ((67|) in 
order to show that, under certain conditions, namely, low intensity (fi <C v) and Lamb-Dicke (rj <^ 1) 
regimes, eigenstates ([6T|) correspond to the blue side band and the eigenstates (|67l) correspond to the 
red side band. Under these conditions is easy to prove that 

m^^V^\l)\e) + \0)\g), (70) 

which is an eigenstate of the operator of the form cr+a^ + (T^a. Thus, when an up ion internal transition 
takes place {\g) — > |e)) the vibrational motion acquires an extra phonon (|0) — > |1)). In general, we 
will have in this approximation 

m 

IV'm) wiVl + m77^d,„|m+l)|e) +^d„|n)|g) (71) 

ri=0 

showing that internal up transitions are accompanied by vibrational up transitions. 

In the case of the eigenstates (|67)) . we have the conditions <C j^, 77 <C 1 and 6 — ly. Under these 

conditions we can approximate the eigenstates as 

m 

IV'm) ~ ^dn\n)\e) - Wl + mri—d,n\m)\g), (72) 

71=0 

thus transitions to the lower internal ion states are associated with increase of vibrational quanta. 

3.3 The blue side band and the red side band by means of the intensity 

In the resonant case, processes that correspond to the blue side band and the red side band can also 
be obtained by means of the laser intensity. If we take (5 = 0, and consider the Lamb-Dicke limit, the 
displacement operator can be expanded in Taylor series and written as 

£)(a)=e""^"-"*'^«l-Haa^-a*a, (73) 
such that the Hamiltonian (|^5|) reads 

H = vn + rifTa; — jycTy (a + a^) , (74) 
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where we have used that ax = + cr_ and (Ty = —i{a+ — cr_). By making now a rotation around the 
Y axis (by means of the transformation exp(i^CTj,)), and going to the interaction picture, we get the 
transformed Hamihonian 



Hi = irjQ 



(75) 



If we take now ly — and we use the rotating wave approximation, the Hamihonian reduces to 

Hi ~ —irfl ( ha- — ct+o^) (76) 



which clearly gives us the blue side band. 

To get the red side band, we must take v — 217 and the Hamihonian we get is 



Hi 



irjfl (c 



(77) 



that clearly implies that when the ion goes from the ground internal state to the internal excited state, 
the vibrational motion losses one phonon and viceversa; i.e., the red sideband. 



4 Solution in different regimes: dispersive Hamiltonians 

The ion-laser interaction may be easily solved in the low intensity regime [351 IMl HOI Ell 112], but 
besides the condition that the laser intensity is much lower than the vibrational frequency, we set 
the condition that the detuning between the laser and the atomic transition frequency is an integer 
multiple of the vibrational frequency. Then some questions arise: Is it possible not to consider integer 
multiples of the vibrational frequency?. Is it possible to solve for high and middle intensities? 

Indeed, we have shown in Section 3 that it is possible to find solutions for any set of parameters; 
i.e., in all regimes j43l; however, the solutions are not general because the set of eigenstates that may 
be found can not expand all possible (general) states. 

It has been shown already that for low intensities it is possible also to consider the ion micromotion 
|26| , and by using Ermakov-Lewis invariant methods ' 301 131[ 132) it was possible to linearize the ion- 
laser Hamihonian when the micromotion was included |44| . Here, we follow Ziiniga et al. [45] and 
show how it is possible to solve the ion-laser interaction in different regimes, including high intensity 
and medium intensity. 



4.1 Different regimes 

In Section 3.1, we showed that the ion-laser Hamihonian (j45p can be casted in the form given by 
expression (j48p by means of the similarity transformation (j46p . Therefore, we have linearized the ion- 
laser interaction in an exact way, by means of a unitary transformation; i.e., both Hamiltonians, i?ion 
and Hion are equivalent. In the following, we will neglect the term because it only represents a 
constant shift of all the eigenenergies. Of course, transformation has to be applied to an initial 
condition for the internal state of the ion and its vibrational motion wave function. Let us assume 
that we have the initial state 

\m) = IWIe), (78) 

where \ia) is a coherent state, and for simplicity we take a a real number (to avoid extra phases 
later, but the calculation may be done for complex a). Then, we have that the initial wave function 
associated with the transformed Hamihonian is 

\m) = T\m), (79) 

where the transformation T is given in (j46p . If we write the initial wave function in terms of 2 x 2 
matrices, we obtain 

D{zrj/2) \ f \ia) \ ^ J_ f |z(a-7y/2)) \ 



\m) 
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Thus, we have changed the comphcated Hamiltonian (j45|) by the Hnear Hamiltonian (|48|) via a unitary 
transformation. The smah prize we have to pay, is that in the initial wave function the coherent state 
is displaced and the ion is initially (in the new frame) in a superposition of ground and excited states. 



4.2 Medium intensity regime (MIR) 

We already considered this case, when we analyzed in Section 3.3 the blue and red side band by means 
of the intensity. In this case, the vibrational frequency is of the order of twice the field intensity 
(Rabi frequency). We also consider the Lamb-Dicke regime; i.e., 77 ^ 1. For simplicity, we will set 
5 = to show the different possibilities we have now. However, it is not difficult to produce effective 
Hamiltonians also in the off-resonance case. In this case, the Hamiltonian (|48p may be casted into (fTTt 
which is a Hamiltonian that has been extensively studied [HIES]; therefore, we will not add more here, 
except the fact that for the medium intensity regime the Hamiltonian (|45p may be exactly expressed 
as a JCM Hamiltonian via a unitary transformation and the rotating wave approximation, without 
extra approximations. 



4.3 Low and high intensity regimes 

In the case of the low intensity regime (LIR) the solution has been known already for several years 
[ini[I3i however, we will treat it with some detail in Section 5. Here, we will show a different method 
that is also valid for the high intensity regime (HIR). Just for the matter of qualitative analysis, let us 
take ^ = 0. Consider now ft i' (LIR) or 51 ^ 1/ (HIR) in equation p8)) . As this Hamiltonian for (5 = 
is equivalent to the atom-field interaction, we can borrow knowledge from such interaction: we know 
that when the field and atomic transition frequencies are very different (in our case, it is translated 
in the equation \i/ — 217| <C t]1'/2, that may happen in either of both regimes, HIR or LIR), atom and 
field stop to exchange energy and we obtain a dispersive Hamiltonian |46| . The same happens in the 
ion-laser interaction, and via an small rotation approach |47| , we will be able to cast Hamiltonian (I48p 
as an effective (dispersive) Hamiltonian. 

By transforming the Hamiltonian (|48|) with the unitary operators 



I.e. 



with ^1,^2 ^ 1, using 



and remaining up to first order in the expansion e^^Be — B + £\A,B] + ^[A. [A,B\\ + ... w 
B + ^[A,B] i47j, we get the effective Hamiltonian 

15 K 

that ior 5 — is known as the dispersive Hamiltonian. Note that, just as in the atom-field case, there is 
no need to transform the (already transformed) initial state ()80p as an small rotation has been applied. 
We can see that in fact ^i, ^2 <C 1 either in the LIR (in this case we have also to consider ry ^ 1) or in 
the HIR (no constrain for 77), which justifies completely the approximation for the above Hamiltonian. 
For the resonant case, S — 0, it becomes diagonal and we can solve it in an easy way. 
In Figure 5, we show a plot for the probability to find the ion in its excited state 

Pe{t) = {^P{0)\fUxp{itn,s)f\e){e\f^exp{-im,s)f\^{0)) (85) 

as a function of time for k = 0, ly = 1, = 0.2 and 77 = 0.1. The vibrational motion of the ion is 
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Figure 5: Plot of Pe{t) as a function of t for k ^ Q, ly = 1, ~ 0.2 and 77 = 0.1. The vibrational motion 
of the ion is considered to be in a coherent state, |ap — 4 and the ion in its excited state. Solid line 
represents the numerical (exact) solution, dashed line the solution from Section 3 and the dot-dashed 
line the solution for the dispersive Hamiltonian. 



— 4, and the internal state of the ion is excited. The three curves 
in the figure correspond to the exact numerical solution (solid line), the solution from Hamiltonian (|45p 
(dashed line) and the solution for the dispersive Hamiltonian ((M)) . We can see excellent agreement 
among the three plots for the LIR. Now, for the HIR in Figure 6, we show a plot also of Pe{t) as a 
function of time for the exact numerical solution (solid line) and our solution from this section (dashed 
line), but now with the parameters for fc = 0, f2 = 1, = 0.2, and rj = 0.1. Again, it may be noticed the 
excellent agreement between both curves. We should stress that there is no other analytical solution 
to compare with, as ours is the first analytical solution in this regime (also in the medium intensity 
regime) . 

The new interaction constants in the effective Hamiltonian (1841) have the forms 



K = —TTT^ -. (86) 



In the resonant case and high intensity regime, i7 3> z^, it is easy to show that 

while in the low intensity regime, fl <^ v, we will have the same Hamiltonian, but x will change to 

1 



Xion Xlow = -2?? n- w -2ry fl. 

If in Equation (|84|) we take the detuning 8 different from zero, we could get the usual blue and red 
side-bands interactions (see for instance Ref. [42 ). This is done by choosing the value (5 = ±i/. The 
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Figure 6: Plot of Pe{t) as a function of t for k ^ Q, ~ 1, i' ~ 0.2 and 77 = 0.1. The vibrational motion 
of the ion is considered to be in a coherent state, jap — 4 and the ion in its excited state. Solid line 
represents the numerical (exact) solution, dashed line the solution for the dispersive Hamiltonian. 



only case in which we can obtain such regimes is the low intensity case, where one can perform the 
rotating wave approximation to the Hamiltonian (j84p . which agrees with the usual procedure for ob- 
taining such blue and red side-band regimes. The high intensity case, fi 3> i^, does not allow such 
side-bands because in the Hamiltonian (j84p the interaction constants multiplying the different terms 
may be of the same order. 



5 Low intensity regime 

If we consider that the Hamiltonian (|38|) corresponds to the wave function the Schrodinger 

equation can be written as 

^§-^\0=H\0- (89) 

Let us examine the transformation to a rotating frame of frequency w, by means of the unitary 
transformation 

T(t) =exp(i|a,t). (90) 
Applying the transformation T, the wave function |^(t)) transforms in the wave function |(/)(i)); i.e., 

mm)) = \m, m 
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and the Hamiltonian transforms in 



H^ = i^^THt)+T{t)HT\t). 



(92) 



Writing the position operator x in terms of the ladder operators, expressions ([39]) and (|40l) . using 
the Baker- Hausdorff formula [261 HSjj and the commutators of the Pauli matrices (|4T|) . the explicit 
transformed Hamiltonian is 



k 



XEo 



g«)(at+a) 



H.C. 



where rj is the Lamb-Dicke parameter. The quantity 

kv = UJ21 — uj 



(93) 



(94) 



is the detuning between the plane wave frequency and the transition frequency of the ion; in other 
words, we are considering that the detuning is a multiple integer of the vibrational frequency of the 
ion. 

We need now to factorize the exponential in the Hamiltonian ([93| . As [a, [a, a^]] = 0, and [a^ , [a, a^]] = 
0, we can use the Baker-Hausdorff formula, and write 



(95) 



Expanding in Taylor series the exponentials that contains the operators, and substituting in the Hamil- 
tonian (j93p . we obtain 



Ht = i^n + jz-cr^ -I- XEoe '''/^ 



E 



(-^?7)"^ 



(a^)"a™-hH.C. 



7 'm ' 



n.m—0 

We go now to the interaction picture, using the transformation 

Tfi-oe = exp [it{i^n + -^<^z)]- 
We apply this transformation, using the following two commutators, 

df 



and 



and we obtain 



E 



[a,/(a,a^) 
'a\f{a,a))] = - 
(-^77)"+™ 



Sat 

dl 

da ' 



7, ' rr?, ' 



(at)"a™e"*(""™"'')''* + H.C. 



Ti.m— 



(96) 

(97) 

(98) 
(99) 

(100) 



where $7 — XEq is the exchange energy frequency of the internal and vibrational states, called Rabi 
frequency. 

The interaction Hamiltonian (llOOp has a diversity of contributions, each contribution oscillates with 
a frequency that is a multiple integer of v. We apply now the rotating wave approximation; as the 
Schrodinger equation is a first order differential equation in time, we have to integrate it once with 
respect to time; this integration brings the sum and the difference of the frequencies to the denominator. 
The terms changing slowly will dominate over the terms changing very fast; so the contribution to 
the Hamiltonian of those fast terms is neglected, and only the slowly changing terms are kept. In this 
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case, the terms that do not rotate quickly are those whose exponent satisfies the relation n — m = k, 
and as we already explained, are those terms that we will keep. This approximation is valid for 



(101) 



As Vl is proportional to the amplitude of the laser electric field, from (|10ip it is clear that this approx- 
imation is valid for low intensity. We have then. 



{-irf) 



2m+k 



(to + k)lml 



Using now the fact that the number states is a complete set 

oo 

where / is the identity operator, we can write 

oo 

(atra'" = (atra"^|j)(j| 



(a^)'=(a^)'"a'" + H.C. 



n=0 



J=0 



(h — to)! ^ 



(n — to)! ' 



that substituted in the Hamiltonian (|102p . gives us 



— ' (to + fc)!TO! (h — to)! 

m=0 ^ ' ' ^ I 



Using the explicit expression for the associated Laguerre polynomials [501 151) 

{n + a)\ 



H.C. 



(102) 



(103) 



(104) 



(105) 



4"n^) = E(-i) 



(n — i)!(Q! + i)! i\ 



we can write 



n! 



Lf (ry2)+H.C. 



(106) 



{h + fc)!^ 

We will consider now processes where only one phonon is exchanged; that means that we must 
take fc = 1 in the Hamiltonian (I106P . We will consider also that the oscillation amplitude of the ion is 
much smaller than the laser frequency; that is, 77 ^ 1; or in other words, we suppose the Lamb-Dicke 
regime. With these two considerations, the Hamiltonian (|106p reduces to (the subscript jc stands for 
Jaynes-Cummings) 

i/jc = ~iT]VL{a^a^ - a(T+). (107) 

In agreement with the considerations made above, the Hamiltonian p07p describes emission and 
absorption of one vibrational excitation, when the atom makes electronic transitions. The first term 
represents the absorption of a vibrational excitation and the transition of the ion from the excited state 
to the ground state. The second term represents the inverse process; the ion goes from the ground 
state to the excited state, annihilating one phonon, and the vibrational state decays in one quanta. 
All this can be clearly seen, if we apply the Hamiltonian (I107P to the correct states. In the first case, 
we have to apply the Hamiltonian to the state 171)16), which represents n vibrational quanta and the 
ion in the excited state |e); we get, 

i7je|n)|e) « |n+l)|5); (108) 

which is the state with n + 1 vibrational quanta, and the ion in the ground state. In the second case, 
the state is given by |n + l)|.g), and when we apply the Hamiltonian we obtain. 



i7je|n+l)|5) cx |7i)|e); 



(109) 
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that is the state with n vibrational quanta, and the ion in the excited state. 

We can repeat all the above procedure now when the laser frequency is greater than that of the 
transition 

kiy ^ uj — UJ21, (110) 
and obtain the Hamiltonian (clearly now, the subscript ajc stands for anti-Jaynes-Cummings) 

Hajc = —irifl{aa^ ~ a^(T+). (HI) 

In the first term, we have the annihilation of one quanta from the vibrational motion and the ion 
internal transition from the excited state to the ground state. In the second term, we have the 
creation of one vibrational quanta and the internal excitation of the ion from the ground state to the 
excited state. In Figure 7, we explain why this Hamiltonian is anti-Jaynes-Cummings type. 
Applying the anti-Jaynes-Cummings Hamiltonian to the adequate states, the previous comments can 





Figure 7: a) A Jaynes-Cummings Hamiltonian implies the ascent (descent) of one ion vibrational 
quantum, and at the same time, the transition from an excited (ground) internal state to the ground 
(excited) state, b) An anti-Jaynes-Cummings Hamiltonian annihilate (creates) one quantum from 
the vibrational motion and transfers the ion internally from the excited (ground) state to the ground 
(excited) state. 

be easily understood. If we apply the Hamiltonian (|llip to the state |n)|e), we get 

i7ajc|n)|e) = |n-l)|g), (112) 

and if we apply it to the state \n — l)!^), we get 

ilajc|n-l)|.g) = |7i)|e). (113) 
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/^From the point of view of the trapped ion, all this means that we can take it to its lowest energy 
vibrational state, alternating successively, and as many times as necessary, the detuning between the 
frequency of the plane wave and the internal frequency of the ion. Again, we can illustrate all this by 
applying the correct Hamiltonian to the adequate state. For that let us consider a vibrational state 
\n) and the ground internal state; if we apply the Hamiltonian ()107p . we get 

H^,\n)\g) = \n-l)\e). (114) 

We apply now the Hamiltonian (|111|) . obtaining 

iJajc|«-l)|e) = |n-2)|5), (115) 

and the ion has lost two quanta of vibrational energy. Repeating successively this procedure, we can 
arrive to the state |0)|5). 

5.1 Adding vibrational quanta 

We now show how to generate nonclassical vibrational states in the low intensity regime. From Equa- 
tion (|102p with fc = 2, we note that if the Lamb-Dicke parameter is much less than one, 77 ^ 1, we can 
remain to the lowest order in the sum, such that we obtain the so-called two-phonon Hamiltonian 

Hi = e[(T^a^ +(j+a^% (116) 

with e — —Q/2. For the study of the dynamics of interest, we need the time evolution described by 
the Hamiltonian ()116|) . The advantage of the interactions of Jaynes-Cummings type consists in the 
fact that the Hamiltonian can easily be diagonalized, and using the same procedure that was already 
used, it is possible to show that 



m— I 

00 



cos ) (In + 2)\g){n + 2\{g\ + \n)\e){n\{e\) 

-ism I ^n,,t] (|7i + 2)|5)(7i|(e| + |7i)|e)(7i + 2|(g|) 



(117) 



^2 

The quantity r2„ is the two-phonon Rabi frequency which is given by 

2ey/{n + l){n + 2). (118) 

Using these results, the time evolution of the quantum state in the interaction picture is easily derived 
for arbitrarily chosen initial conditions. We have 

\^{t))^Uj{t)\^{0)). (119) 

If we consider as initial state the ion in its excited state |e) and the vibrational state a coherent state, 
we can find the atomic inversion, (that we recall that it is defined as the probability to find the ion in 
its excited state minus the probability to find it in the ground state). Using (|117p and (|119p . we get 



W(t) = exp (-lap) V cos(2nnt). (120) 

Tl = 



2n 



We plot this function in Figure 8 as a function of the scaled time t — et. The interaction gives rise to 
a quasi-regular evolution of the atomic inversion, unlike the case of one phonon resonance. This can 
be used for several purposes, among them, to add excitations to the vibrational state. In Figure 8, it 
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Figure 8: Plot of the atomic inversion, W{t) = P2 — Pi, the probabiUty to find the ion in its excited 
state minus the probability to find it in the ground state, for an ion initially in its excited state and 
the vibrational state in a coherent state, with a — 5. 

can be seen that if initially the ion is in its excited state |e), after an interaction time t — ir/e the ion 
ends up in its ground state \g), giving all its energy to the vibrational state by adding two vibrational 
quanta. Moreover, the effect of having the ion in its excited state and after an interaction time having 
it in the ground state is shared by all vibrational states, not only when it is prepared in a coherent 
state. To illustrate this fact, we show in Figure 9 the atomic inversion for a thermal distribution 

n=0 ^ ' 

with n the average number of thermal phonons. 
Consider again the initial state 

l*(0)) = |a)|e). (122) 

Combining Equations (jll7p and (|122l) . and using a compact operator representation of the Jaynes- 
Cummings dynamics, we arrive at 




(123) 



In order to derive illustrative analytical results, in the following we will apply the approximation 

y^a2(^t)2«^+| (124) 

although this approximation represents a Taylor-series expansion for large eigenvalues n of the operator 
n, the error is already small for small n-values. For example, for n = 1 the relative error is only 0.02. 
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Figure 9: Plot of the atomic inversion, W{t) = P2 — Pi, the probabiUty to find the ion in its excited 
state minus the probabihty to find it in the ground state, for an ion initially in its excited state and 
the vibrational state in a thermal distribution with n — 2. 

Based on this approximation, one may simplify Equation (jl23p as 

l^'(i)) « cos[Ai(n + 3/2)] |a) |2) - ^(y^)^ sin[A<(7i + 3/2)] (125) 
Choosing a particular interaction time t = r, according to 

T = tt/A, (126) 

we obtain for the vibrational state vector 

\¥^\r))^^^iV^r\-a), (127) 

where we have introduced the subscript "w" (vibrational) to note that we are not taking into account 
anymore the state \g). Moreover, the subscript "+" and the superscript "(fc)" are used to indicate the 
process of adding (two) vibrational quanta and the number of such interactions respectively. 
The quantum state (jl27p may serve as the initial state for a second interaction with an ion that is 
prepared in the same manner as the first one. For the same interaction time, t — t, after the second 
interaction (that is completed at time 2t), the vibrational state is 

\^^^\r)), = -iV^)^\a). (128) 

By repeating the process k times, one finally obtains for the quantum state, at the time tk, after 
completing k interactions, 

(r)), -^^■(V'^)''l(-l)'«)■ (129) 

After many interactions, the state |\I'^^(r))i, exhibits a strong sub-Poissonian character, because while 
one is adding two excitations per interaction, at the same time one is keeping the width of the distri- 
bution constant. 
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(k) 

The excitation distribution Pn , after k interactions, is easily found to be related to the number 
statistics Pi"'' of the initial state la) via 



p(fc) _ p(0) 



(130) 



This result clearly shows that the number statistics is only shifted but retains its form. 

One way of studying the properties of the states being generated is through the Mandel Q-parameter 

[5^ . which is defined by 



(131) 



and where 



if g 



> 0, super Poissonian distribution 

= 0, Poissonian distribution (coherent state) 

< 0, sub-Poissonian 

= — 1, number state. 



In this case the Mandel Q-parameter is given by 



(132) 



2k 



1, 



(133) 



and as the number of interactions increases, the Mandel Q-parameter approaches the value —1; i.e., 
the state acquires maximum sub-Poissonian character. The sub-Poissonian effect of the vibrational 
wave function becomes more significant with increasing number of interactions. 



5.2 Subtracting vibrational quanta 

It is straightforward to show that, opposite to the case in which we add two phonons per interaction 
by initially having the ion in its excited state, two phonons may be removed per interaction by initially 
having the ion in its ground state. We therefore expect that instead of squeezing the phonon number 
distribution (giving a sub-Poissonian character to the distribution), in the case in which we subtract 
phonons, it should be broadened; however, the conjugate variable to the number operator; i.e., the 
phase operator [49] should have less fluctuations; a squeezing of the phase distribution then will happen. 
It is not difficult to show that the phonon distribution when we subtract two phonons per interaction 
will be given by ^12] 

Pi'' - PiX, (134) 
from which we can calculate the phase distribution plotted in Figure 10 for several number of interac- 
tions. 



5.3 Filtering specific superpositions of number states 

If instead of one laser, we assume two lasers driving the ion, the first tuned to the jth lower sideband 
and the second tuned to the mth lower sideband, we may write E^^'>{x,t) as 

E(-\x,t) = Eje-'(''^^-'^^'+^'''>* + E„,e-'^''"^^-'^^'+""''>\ (135) 

where if to = 0, it would correspond to the driving field being on resonance with the electronic 
transition. The position operator x may be written as before, 

ksX = ris{a + a^), (136) 

where ks, s = j,m are the wave vectors of the driving fields and 



V(0|A£2|0) 

fjs = 2Tr^ (137) 
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p('^)(o) 




Figure 10: Subtracted phase distribution for an initial coherent state a = 45 and (a) k ~ 0, (b) 
k = 900, (c) k = 1000 and (d) k = 1012. 



are the Lamb-Dicke parameters with s = j, m. 

In the resolved sideband limit, the vibrational frequency v is much larger than other characteristic 
frequencies, and the interaction of the ion with the two lasers can be treated separately using a 
nonlinear Hamiltonian [16] . The Hamiltonian (|106p in the interaction picture can then be written as 



Hi = 



H.C. 



(138) 



where Lf^\rj'^) are the operator- valued associated Laguerre polynomials, the Q's are the Rabi frequen- 
cies and h — a) a. The master equation which describes this system can be written as 



dp r 

— = -i[Hi, p] + 77 ( 2o-+pcr_ - a^p - pa, 



dt ' 2 

where the last term describes spontaneous emission with energy relaxation rate F, and 



i 1 
P=2 



(139) 



(140) 



accounts for changes of the vibrational energy because of spontaneous emission. Here rjE is the Famb- 
Dicke parameter corresponding to the field ()135p and W{s) is the angular distribution of spontaneous 
emission |16j . 

The steady-state solution to Equation ()139p is obtained by setting dp/dt ^ 0, and may be written as 



Ps = \ilJs)\9){ips\{9\, 



(141) 
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where \g) is the electronic ground state and {ips) is the vibrational steady-state of the ion, given by 



l^s) = 0. 



(142) 



For simplicity, we will concentrate in the j = 1 and m = case (single number state spacing) for which 
Equation (|142p is written as 



(143) 



Note that Hi\l)\ips) — 0, so that ion and laser have stopped to interact, which occurs when the ion 
stops to fluoresce. For the j — I and fc = case, and assuming L''j^\ril) ^ and Lk{r]Q) ^ for all 
k, one generates nonlinear coherent states [13]. However, by setting a value to one of the Lanib-Dicke 
parameters such that, for instance, 

LM) = 0, (144) 
for some integer q, we obtain that, by writing \^jjs) in the number state representation. 



(145) 



Tl = 



(the argument of ipg denotes the condition we apply; i.e., in Equation (|145p . the condition is on 770) 
where 



("0^/^ n 



'.LIMY 



(0) 



= 1, 



and 



is the normalization constant. 
If instead of condition (|144p . we choose 



we obtain the wave function 



n— p+1 



where now 



floe-"'/' 



n—p— 1 



^(1) 



= 1, 



and 



n— p-f 1 



(146) 
(147) 

(148) 
(149) 



(150) 
(151) 
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Combining both conditions, (|144|) and (|148|) . one would obtain for q > p, 



IMvo.m))^^ J2 C^n^n), (152) 

n— p+1 



with 



N'o^^ E (153) 

n— p+l 

In this way, by setting the conditions (I144p . (|148p or both, we can engineer states in the following 
three zones of the Hilbert space: (a) from |0) to \q), (b) from \p+ 1) to |oo), or (c) from \p+ 1) to \q). 
In the later case, by setting q = p + I, generation of the number state \q) is achieved. 
We should remark that by selecting further apart sidebands one would obtain a different spacing in 
Equations p45p . (|149l) and (|152p . For instance, by choosing j = 2 and fc = one would obtain only 
even or odd number states in those equations (depending in this case on initial conditions, and W{s), 
the angular distribution of spontaneous emission). Also, it should be noticed that one can use the 
parameters j — m + 1 and k = m (with to ^ 0) (in the single number state spacing case) to extend 
the possibilities of choosing Lamb-Dicke parameters. Lamb-Dicke parameters of the order of one (or 
less) are needed (for conditions (|144p and (|148l) ). which can be achieved by varying the geometry of 
the lasers. For example, by setting rjo = 1, we have Li(j]q = 1) = 0, and therefore we obtain the qubit 

\Mvo = 1)) = , ^ , (|0) - ^^^Il>) , (154) 

where by changing the Rabi frequencies, one has control of the amplitudes. Finally, note that we could 
have also chosen to drive the gth upper sideband instead of the kth lower sideband in Equation (|135p 
with basically the same results. 



5.4 A^OOA^ states 

Nonclassical states have attracted a great deal of attention in recent years, among them (a) macroscopic 
quantum superpositions of quasiclassical coherent states with different mean phases or amplitudes 
[SSHM]) (b) squeezed states [551 [21], (c) the particularly important limit of extreme squeezing; i.e., 
Fock or number states, and more recently, (d) nonclassical states of combined photon pairs also called 
NOON states [56l [57] . It is well known that these multiphoton entangled states can be used to obtain 
high-precision phase measurements, becoming more and more advantageous as the number of photons 
grows. Many applications in quantum imaging, quantum information and quantum metrology [58j 
depend on the availability of entangled photon pairs because entanglement is a distinctive feature of 
quantum mechanics that lies at the core of many new applications. These maximally path-entangled 
multiphoton states may be written in the form 

INOON)^, = i= |0), + |0), |iV),) . (155) 

It has been pointed out that A^OOA^ states manifest unique coherence properties by showing that they 
exhibit a periodic transition between spatially bunched and antibunched states when undergo Bloch 
oscillations. The period of the bunching/antibunching oscillation is N times faster than the period of 
the oscillation of the photon density [55^ . 

The greatest N for which NOON states have been produced is iV = 5 [56]. Most schemes to generate 
this class of states are either for optical [SHI [57] or microwave [50] fields. In this contribution, we would 
like to analyze the possibility to generate them in ions [SHI [Ml [HI [121 [13]; i-e., NOON states of their 
vibrational motion. We will show that they may be generated with N = 8. 
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Figure 11: Possible situations we can have if we filter number states with the proposals of this section. 



5.4.1 Ion vibrating in two dimensions 

We consider an ion in a two-dimensional Paul trap |61) , and we assume that the ion is driven by a 
plane wave 

E^-\x,y,t) ^ Eoe-"'-''^'^+''yy+'^^\ (156) 
with kj^j = x,y the wavevectors of the driving field. The Hamiltonian has the form 

H ^i^^a^ax + Vyalay + -— CTz 

(157) 

|e-4''-("-+'ii)+Ha+ + H.C.} + VLy |e-''[''«(^«+4)+'^*]<7+ + H.C.} . 
where we have defined the Lamb-Dicke parameters 



,..2,vMC. „,.2,vSWK, ,158) 

and redefined the ladder operators according to 

Kx ^ Tj^ia^ + al), k i^ndV ^ Vyi^y + al). (159) 

In the resolved sideband limit, the vibrational frequencies i^x and Uy are much larger than other 
characteristic frequencies and the interaction of the ion with the two lasers can be treated separately 
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using a nonlinear Hamiltonian [T6l [17] . 

We consider that the ion is trapped in the x - axis; i.e., Q.x ^ ^ and U,y = 0; then 

= v.ala, + + f!, |e-4''^(a.+a^)+"*]a+ + H.C.} . (160) 

We write W21 = w + (5, where (5 is the detuning, to obtain 

^^^^<7, + n., |e-4"^(a^+-l)+Ha+ + H.C.} . (161) 



We transform the Hamiltonian to a frame rotating at uj frequency by means of the transformation 



and we get 



(162) 



(163) 



Using the Baker-HausdorfF formula [48] . and expanding the exponentials in Taylor series, we cast the 
Hamiltonian to 



Hx = i^xal&x + -dz + 



e 2 



Going now to the interaction picture. 



e 2 



E 



aT"a;"a+ + H.C. 



(164) 



(165) 



We consider now the low- intensity regime; i.e., Qx << Vx, and we apply the rotating wave approxima- 
tion, to get 



Hjx — ^x 



2 00 , -.-zn 

n! (fe -|- nj! 



by substituting aj."a" = , multiplying by , and rearranging terms 



e 2 



(166) 



(167) 



where we have identified L| (77^) = ^ ^ ^^J''"'"^ (n-Hfc)T(ft-n)! ' ^^^'^ ^'^^ associated Laguerre polyno 

n— 

mials; so that finally 

Hix = nx {f^ (n) a^a+ + ai>^ {h) a_) , 



where 



fx (n) - e "^"^ i-iVx) 



(n + Zc)! 

We write the Hamiltonian, given by expression (jl68p . in the following matrix form 



Hix 



nxfi^in)a'x 



which, by using the Susskind-Glogower phase operator i62| . we can be written as 

1 



Hix - 



1 

yjfe 



Hi 



v: 



(168) 
(169) 

(170) 
(171) 
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where we have introduced and defined 



n^f^ in) ^fi^ 



(172) 



The evolution operator for this last transformed Hamiltonian, Uix = e ^Hi^t ^ calculate easily. 

For this we need 



H 



2m 



2m 



2ni 



(173) 



and 



jj2m+l ^ 



^2m+l\ 





\S^{n)\ 



2m+l 



Then 



and therefore 



/ ■,\2m 



.,x2m+l 



( TT2m+l 

(2m +1)! ' 



/ 

(174) 
(175) 



Uix (i) 



cos (n^t I (n) I -i {-if sin (n^t \ (n) | V? 



V 



iV^^ {if sin n,t I {n) I J dial'' cos O.t I (n) I J ald^f 



Finally, as \j a\a}x^ = \f^^^^ write 
/ 

Ulx{t) = 



(176) 



cos 



v 



Oqj^ I 

-(z)'=+V/sin 



(n+fc)! 



COS 



f^.i {n) 



(h+k)\ 



yk 



J 



(177) 

Now, wc consider as initial state of the ion a number state \n) for the vibrational motion and the 
excited state |e) for the internal states; i.e.. 



i^(o)) = (';;^ 

The probability, after the time of finding the ion in its internal excited state is then 

oo / 

Pe {t) = Y.{e\{m\i,{t)){m\m)\e) = cos' {n)\ 

m=0 \ 

It is clear that after a time 



{n + fc)! 



77 1 
i0 = 7T 



2n,\f^{n)\y {n + ky. 



(178) 



(179) 



(180) 
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1 2 3 4 5 6 7 

Figure 12: The probability, after the time r, of finding the ion in its internal excited state when initially 
the ion was in its excited state and the ionic vibration was a number state with n = A and rj^ =0.1. 

the probability to find the ion in its internal excited state is 0, so at that time the ion is in its 
internal ground state with probability 1. This situation is obviously repeated periodically; every 
2i + 1, j = 0, 1, 2, ... times to, the ion will be in its ground state. We plot in Figure 12 the probability 
for the case when fc = 4 and the ionic vibration is initially in a number state with n = 4. As can bee 
seen from Hamiltonian (jl68p . when the probability of finding the ion in the excited state goes to zero, 
the ion is giving four phonons to the vibrational motion. Now, if we consider the ion initially in its 
ground state, the probability to find it in the ground state at the same time tg is also zero. In this 
case the ion removes four phonons of the vibrational motion. 

If we consider now that the ion is trapped in the y - axis; i.e., fly and fix — 0, we get exactly the 
same expressions and the same results with the variable y instead of x. 

5.4.2 Generation of NOON states 

By starting with the ion in the excited state and the vibrational state in the vacuum state; i.e., 
|0)x|0)j,, if we set r]y = 0, after the time Tp when the probability to find the ion in its excited state 
is zero (meaning that the ion, by passing from its excited to its ground state, gives 4 phonons to 
the vibrational motion), we can generate the state |4)3,|0)j^. Repeating this procedure (with the ion 
reset again to the excited state, via a rotation), but now with rj^ — 0, four phonons are added to the 
2/- vibrational motion, generating the two-dimensional state |4)2,|4)j^. 

Therefore, if we consider the ion initially in a superposition of ground and excited states, and the 
|4)a;|4)j, vibrational state; i.e., 

I^^t) = -^(|e) + |5))|4),|4)„ (181) 
for rjy — and tg, the state generated is 

|V'r,„=o) = ^(|e)|0). + |<?)|8).)|4),. (182) 
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Now, we consider this state as initial state for the next interaction with rj^ = and still the interaction 
time Tp, to produce 

\^,.=o) = -^(|e)|0),|8), + |g)|8),|0),). (183) 
Next, the ion is rotated via a classical field (an on-resonance interaction) such that the state 

\^r) = -I [|e)(|0).|8), ~ |8),|0),) + |5)(|0).|8), + |8),|0),] (184) 
is obtained. Finally by measuring the ion in its excited state we produce the iVOOTV state 



liVOOTV,) = -^(|0),|8), - |8),|0),), (185) 



and if the ion is measure in the ground state, also a NOON state is produced: 



|iV00Ar^) = i=(|0),|8), + |8),|0),), (186) 



5.5 Measuring squeezing 



In this section, we have shown how to generate nonclassical states of the vibrational motion of the 
ion. The question arises: Can we measure nonclassical features in this interaction? Here we propose a 
method to measure squeezing. In order to achieve this we need to be able to measure quantities like 

(1) -(a)+c.c., {X^)^{a^) + {[a^]^)+2{h) + l. (187) 

Below we will show how it is possible to measure such quantities by utilizing an atom as a measuring 
device. Consider the Hamiltonian (109), which in matrix form may be written as 

H = -if]VL{a^&- - a+a) = -iryl] (^^^ ) ' '''^^^^ 

We can re-write Hamiltonian (|188p with the help of Susskind-Glogower operators [62] as 



where 

R=i\ I ) (190) 

as 



^ ^ 

with V — ,} , , a. Note that R^R = 1, but RR^ ^ 1. This allows us to write the evolution operator 



U{t) = R ( 'JJ*^{5±I) - ^){^) ) i?t (191) 

We are neglecting a term |0)(0| in the above evolution operator (in the element "2,2"), that however 
will not affect the measurement of squeezing as we will consider the ion in the excited state. We 
consider the vibrational wave function in an unknown state, such that the initial state of the system 
is 1^(0)) — 16)1-0^(0)); the average of the operator <7+ is given by 

(ct+) = {tl}^{Q)\cos{rinWn + l)V^ sin(r;f7tVn + l)IV't.(O)) (192) 
= \{'^vmy^ (sin[?7mA+(n)] - sin[77r!tA_ (n)]) |V'„(0)), 



33 



where 

A+(7i) ^ \/n + 2 + y/n + l, A_ (h) ^ \/n + 2 - y'n + 1. 
By integrating (|193p by using a Fresnel integral [51] 



J dTTsm{T^/A) sm{BT) ^ — f cos 
such that (with rjflt = T) 



dTT siniry A){a+) = --(^^(0)|l/t(^^ „^,)|^^(o)) 



with 



and 



.(1) ^ AA+{n) UA_ 
4 V 2 



cos 



(n) 



+ sin 



AA2^(n) 



.(1) AA_(n) /^rZ/ 

79 = : \ cos 

'2 4 V 2 I 



"AA^(n)' 


+ sin 


'AAl{n)' 


) 


4 




4 





(193) 
(194) 

(195) 
(196) 
(197) 



Now we use the approximation [63j |64] \/(n + 2)(n + 1) w n + 3/2 that is valid for large photon 
numbers; we then can write A^(n) « 4fi + 6 and A^(n) « 0. By setting A = An, we obtain 



7^ ^ « (Vn + 2 + + l)7r cos [(4n + 6)7r] = \/27r^A+(n) 



and 



72 ' « V2^2^_(n), 
so that the integral transform (|195p becomes 



drTsin(rV^)(a+) = V27r2(V'^,(0)|T>^\/^TT|V«(0)) 
= V2^2(V'.(0)|a^|7A.(0)). 



(198) 
(199) 

(200) 



To measure ('0F(O)|[a^]^|V''F(O)) it is necessary a two-phonon transition; in this case 
H2 = X'^^^R^ 







v/(n + l)(n + 2) 
v/(fj + l)(n + 2) 



[i^t]^ 



(201) 



where A^^-* is the interaction constant in the two-phonon case. One can find the evolution operator that 
will be given by an expression similar to ()19ip . just changing y/n + l — >■ ^y(fl+Y)(fl + 2^ , V ^ V"^ 

and [^^]^- It is then easy to calculate the average of ct''^\ which is given by 



-(2)\ 



-i(V'F(0)| cos \\^^h^{h + l){n + 2) \ [V^]^siM^h^{h + l){n + 2) \ |Vf(0)) 
-\{i^Fm[y^f (sin[Ai^+(n)] - sin[AtL(n)]) |Vf(0)) (202) 



with 
and 



5+{n) = V(n + 4)(n + 3) + ^J{h + 2){n+l) w 2n + 5, 



5-{n) = yj{n + 4)(n + 3) - [h + 2){n + 1) « 2. 



(203) 
(204) 



34 



Again by (Fresnel) integration of the above expression 

poo 

/ dTT sin(rVA)(af) = -»47r2(^^(0)|[yt]2(^(2) _ ^(2))|^^(0)) 
Jo 

with 



(2) _ Ad+{h) ttA 



and 



7i 



.(2) _ AS^ jh) hKA 
'^2 - 4 V 2 



cos 



A6l{h) 



AS^h) 



+ sm 



ASlifi) 



Adl{h) 



By choosing the value A = Svr, we obtain 



(205) 



(206) 



(207) 



dTTsm{T^/8TT){a^_^^) 



-iSn^ {iPf (0) \[V^f ^(h + l){h+\)\^F (0)) 
Therefore, squeezing may be measured via this scheme. 



(208) 



6 Ion- laser interaction in a trap with time-dependent frequency 

In this section, we study the problem of an ion trapped with a frequency that depends on time and 
interacting with a laser beam |65[ 166] . Using unitary transformations, we show that this system is 
equivalent to a system formed by a two levels subsystem with time dependent parameters interacting 
with a quantized field. The procedure to build the Hamiltonian for this case is exactly the same that in 
the time independent frequency case (Section 3), but we have to keep in mind that now the frequency 
is time dependent. 
The Hamiltonian is 



and then the Schrodinger equation can be written as 

^^^m) = H\m)■ 

To solve the problem, we make the transformation 
where 

?SD {t) = exp ■ 



zln [p{t)y/iAi\ {xp + px) 







1 exp 


[ 2p(<) J 



and we have also the Ermakov equation 



(209) 

(210) 

(211) 
(212) 

(213) 



as an auxiliary equation. We apply the unitary transformation (|212[) to the Hamiltonian (|209l) : i.e. 
we must calculate the expression 



ot 



(214) 
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We find, 



+ n {exp [-1 {kp{t)y/v^x^ - Lot)] a- + H.C.} 



(215) 



with r2 — XEo, the Rabi frequency. 

Using the Ermakov invariant, the time dependence of the trap has been factorized; the time depen- 
dence is implicit in p(t). We go now to a frame rotating at frequency w, by means of the unitary 
transformation 

(216) 



T^{t) = exp ( -ujta, 



The Hamihonian is transformed to 



(p^ + t'ox^) + 7:i^2i - io)az 



1 



+ n{t) {exp [-i{a + a^)r]{t)\ O- + H.C.} . 



(217) 



Denoting the detuning frequency between the laser and the ion by (5 = l>J2\ — and the characteristic 
frequency of the time dependent harmonic oscillator by 



w(t) 



1 



we get 



Cj{t){h + 1/2) + -cr^ + f](t) {exp [-i(a + a^)r\{t)\ <t_ + H.C.} . 



The time dependent Lamb-Dicke parameter is 
with 

Vo = k\ 

where k is the wave vector of the laser beam. Comparing with the Hamiltonian 
(|219l) is equivalent, but with all the parameters depending on time. 

6.1 Exact linearization of the system 

We call linearization of the system the process to reduce the exponent of the ladder operators a and 
to the first power, without using approximations. To this end, we make the transformation 



(218) 

(219) 

(220) 
(221) 

the Hamiltonian 



where R{t) is given by 



R{t) = exp --((T+ - cr_) 



exp 



_,!|)(a + aty^ 



Transforming the Hamiltonian, we get 

HR = iyon + na, + l^^+i [al3(t) - a^f3*{t)] | {<t+ + (7_) 



(222) 



(223) 



(224) 



36 



with 



(225) 



2 2 

The term ilj{t)/2 has not been considered, because it is only a phase, and when the observable mean 
values are taken, it disappears. 



With the transformation (j223p we have achieved our goal: linearize the Hamiltonian without any type 
of approximation. The rotating wave approximation is not used, and this leaves open the possibility to 
consider different intensity regimes. No assumption has been made about the Lamb-Dicke parameter 
ri[t). It is also valid for any type of detuning and for any time dependence of the frequency of the trap. 
It is also important to remark that we have not imposed any condition in the time dependence of the 
frequency of the trap; in principle, this frequency can assume any temporal form. For a Paul trap, the 
more general form is 



2gcos 2t, 



(226) 



and the Hamiltonian (|224[) is the ion-laser interaction with micromotion included. Also, this Hamilto- 
nian gives us the freedom to consider arbitrary time dependent frequencies. For instance, if we consider 
a sudden change in the trap frequency, we would generate squeezed states for the vibrational wave 
function. 



6.2 Squeezed states by changing the trap's frequency 

If we consider no interaction with a laser; i.e., $7 = 0, the Hamiltonian for the ion with an arbitrary 
(trap) frequency is simply |67] 

H^]^[p'+i^\t)x^] (227) 
and we have seen that the transformation (j212[) produces the Hamiltonian 



2,0 
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Figure 13: Trap's time dependent frequency as a train of steps. 



H = 



(228) 
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which is in an integrable form. If we consider the vibrational motion state to be in a coherent state, 
la), then the evolved wave function reads 



where 



(229) 



(230) 



is the so-called Lewis-Ermakov invariant. If we consider i>{t) to have the form as in Figure 13, we 
can solve numerically the Ermakov equation, and obtain the auxiliary function, p, which we plot in 
Figure 14; note that we will have a series of maximums and minimums. Because the transformation 
Tgj^(t) depends on p, we can analyze from this figure the form in which the vibrational wave function is 
modified. Note from ()212p that the second factor involves the derivative of p, which for the maximums 
and minimums is zero, making this exponential equal to one. Therefore transformation pi2|) may be 
written for such times as 




Q 



50 



100 



Figure 14: Solution for the Ermakov equation for a time dependent frequency as shown in Figure 13. 



TsuitAi) = exp • 



i In [p{tM)y/i^] {xp + px) 



(231) 



i.e., a squeeze operator. Because Tsd{^) — 1 for this form of the time dependent frequency, if initially 
we start with a coherent state, squeezed states will be produced at times Im- 



7 Nonlinear coherent states and their modeUng in photonic 
lattices 

Nonlinear coherent states [551 US] may be generated in ion traps [T3] and in this way realization of a 
quantum-mechanical counterpart of nonlinear optics has been achieved. For an appropriate laser-beam 
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propagation geometry which affects only the dynamics in one vibrational mode of frequency in the 
rotating-wave approximation, the Hamiltonian describing the effect of the Raman laser drive on the 
dynamics of the vibrational mode is given by |13j 

B ^^gk{n){iT]af ^^.G., (232) 

with (n|gi;(ri)|n) = {^Xhy. 

we consider the initial vibrational wave function to be 
in the vacuum state, after application of the evolution operator to this initial condition, a nonlinear 
coherent state is obtained. If we write the solution for the above Hamiltonian in the form 

oo 

IV-W) -^u„(t)|n), (233) 

n=0 

and insert it into the Schrodinger equation, we obtain the following semi-infinite system of differential 
equations for the amplitudes Unit) (for fc = 1) 

and 

= ^{9i{n)VnUn^i + gi{n + l)Vn+ lu„+i), n > (235) 
This system may be modeled in classical optics by light propagation in optical lattices as evanes- 




Figure 15: A semi-infinite waveguide array where fk = gi(k)\/k is the corresponding coupling constant. 

cently coupled waveguides have emerged as a promising candidate for the realization of an ideal, 
one-dimensional lattice with tunable hopping [701 171] . The system of differential equations given in 
(j234p and (|235p may be produced in a setup of waveguide arrays as the one in figure 15, therefore 
modeling the ion laser Hamiltonian (|232p . 

8 Conclusions 

We have studied the ion-laser interaction, producing solutions that do not require too many approxi- 
mations and allowing to reach ranges of parameters, such as Lamb-Dicke, detuning and intensity, not 
reached before. Exact linearization in the case of time dependent frequency has been shown by using 
Ermakov-Lewis invariant methods. In the case in which the frequency is considered constant, we have 
shown that it is possible to find exact eigenstates, access blue and red sidebands with (a) in the low 
intensity regime without considering the rotating wave approximation and (b) in the on resonant case. 
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by adjusting the intensity to half the trap's frequency, i.e., solely by intensity manipulation. We have 
also shown how a special class of nonclassical states, namely, NOON states may be produced in these 
systems and how the ion laser interactions may be mimicked by a particular evanescent coupling of 
waveguides. 
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